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Abstract. We give a systematic discussion of the contributions 
to the cosmological energy density in axions from vacuum re- 
alignment, string decay and wall decay. We call these the cold 
axion populations because their kinetic energy per particle is at 
all times much less than the ambient temperature. In case there 
is no inflation after the Peccei-Quinn phase transition, the value 
of the axion mass for which axions contribute the critical energy 
density for closure is estimated to be of order 6 • 10^^ eV, with 
large uncertainties. It is emphasized that there are two groups 
of cold axions differing in velocity dispersion by a factor of order 
10^. 



1. Introduction 

The axion jl], ||] was postulated approximately twenty years ago to explain 
why the strong interactions conserve the discrete symmetries P and CP 
in spite of the fact that the Standard Model of particle interactions as a 
whole violates those symmetries. It is the quasi-Nambu-Goldstone boson 
associated with the spontaneous breaking of a Upq{1) symmetry which 
Peccei and Quinn postulated. At zero temperature the axion mass is given 
by: 
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where Va is the magnitude of the vacuum expectation value that breaks 
Upq{1) and is a strictly positive integer that describes the color anomaly 
of Upq{1). The combination of parameters fa = Va/N is usually called 
the axion decay constant. Axion models have degenerate vacua [||, 
Searches for the axion in high energy and nuclear physics experiments have 
only produced negative results. By combining the constraints from these 
experiments with those from astrophysics ^, one obtains the bound: 
rua < 10-2 eV. 

The axion owes its mass to non-perturbative QCD effects. In the very 
early universe, at temperatures high compared to the QCD scale, these 
effects are suppressed and the axion mass is negligible. The axion mass 
turns on when the temperature approaches the QCD scale and increases 
till it reaches the value given in Eq.([|) which is valid below the QCD scale. 
There is a critical time ti, defined by ma{ti)ti — 1, when the axion mass 
effectively turns on ||^ . The corresponding temperature Ti ~ 1 GeV. 

The implications of the existence of an axion for the history of the 
early universe may be briefly described as follows. At a temperature of 
order Va, a phase transition occurs in which the U pq(1) symmetry becomes 
spontaneously broken. This is called the PQ phase transition. At that time 
axion strings appear as topological defects. One must distinguish two cases: 
1) inflation occurs with reheat temperature higher than the PQ transition 
temperature (equivalently, for the purposes of this paper, inflation does not 
occur at all) or 2) inflation occurs with reheat temperature less than the 
PQ transition temperature. 

In case 2 the axion field gets homogenized by inflation and the axion 
strings are 'blown away'. When the axion mass turns on at ti, the axion 
field starts to oscillate. The amplitude of this oscillation is determined 
by how far from zero the axion field is when the axion mass turns on. 
The axion field oscillations do not dissipate into other forms of energy 
and hence contribute to the cosmological energy density today Such 
a contribution is called of "vacuum realignment" . Note that the vacuum 
realignment contribution may be accidentally suppressed in case 2 because 
the axion field, which has been homogenized by inflation, may happen to 
lie close to zero. 

In case I the axion strings radiate axions |@, || from the time of the 
PQ transition till ti when the axion mass turns on. At ti each string 
becomes the boundary of N domain walls. If = f , the network of walls 
bounded by strings is unstable |lO) and decays away. If > 1 there is a 
domain wall problem |^ because axion domain walls end up dominating the 
energy density, resulting in a universe very different from the one observed 
today. There is a way to avoid the domain wall problem by introducing 
an interaction which slightly lowers one of the A^ vacua with respect to the 
others. In that case, the lowest vacuum takes over after some time and the 
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domain walls disappear. There is little room in parameter space for that 
to happen and we will not consider this possibility further here. A detailed 
discussion is given in Ref . . Here we assume = 1 . 

In case 1 there are three contributions to the axion cosmological energy 
density. One contribution |Q, | , 12, 13|, |l^ is from axions that were radiated 
by axion strings before ti; let us call it the string decay contribution. A 
second contribution is from axions that were produced in the decay of 
walls bounded by strings after ti |l^, |l^, 11 ; call it the contribution 
from wall decay. A third contribution is from vacuum realignment . To 
convince oneself that there is a vacuum realignment contribution distinct 
from the other two, consider a region of the universe which happens to be 
free of strings and domain walls. In such a region the axion field is generally 
different from zero, even though no strings or walls are present. After ti, the 
axion field oscillates implying a contribution to the energy density which is 
neither from string decay nor wall decay. Since the axion field oscillations 
caused by vacuum realignment, string decay and wall decay are mutually 
incoherent, the three contributions to the energy density should simply be 
added to each other [p^ . 

We will see below that the axions from vacuum realignment, string 
decay and wall decay are cold, i.e. their effective temperature is much 
smaller than the temperature of the ambient photons. In addition to cold 
axions, there is a thermal axion population whose properties are derived by 
the usual application of statistical mechanics to the early, high temperature 
universe jlTt . The thermal axions have an effective temperature of order 
that of the ambient photons. The contribution of thermal axions to the 
cosmological energy density is subdominant if iria 10~^ eV. We concern 
ourselves only with cold axions here. 

The next three sections discuss the contributions to the cosmological 
axion energy density from: 

1. vacuum realignment 

a. "zero momentum" mode 

b. higher momentum modes 

2. axion string decay 

3. axion domain wall decay. 

The basis for distinguishing the contribution from vacuum realignment |^ 
from the other two is that it is present for any quasi-Nambu-Goldstone field 
regardless of the topological structures associated with that field. Contri- 
butions la and lb are distinguished by whether the misalignment of the 
axion field from the CP conserving direction is in modes of wavelength 
larger (la) or shorter (lb) than the horizon size ti at the moment the ax- 
ion mass turns on. Contribution 2 is from the decay of axion strings before 
the QCD phase transition. Contribution 3 is from the decay of axion do- 
main walls bounded by strings after the QCD phase transition. In case 
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1, all contributions are present. In case 2, contributions lb, 2 and 3 are 
absent or at any rate much suppressed; contribution la is present but, as 
was already mentioned, it may be accidentally small. 

In section 5, we estimate the total axion cosmological energy density. 
In section 6, we point out that there are two groups of cold axions (pop. I 
and pop. II) differing in velocity dispersion by a factor of order 10'^. 



2. Vacuum realignment 

Let us track the axion field in case 1 from the PQ phase transition, when 
Upq{1) gets spontaneously broken and the axion appears as a Nambu- 
Goldstone boson, to the QCD phase transition when the axion acquires 
mass. For clarity of presentation we consider in this section a large region 
which happens to be free of axion strings. Although such a region is rare in 
practice it may exist in principle. In it, the contribution to the cosmological 
axion energy density from string and domain wall decay vanishes but that 
from vacuum realignment does not. In more typical regions where strings 
are present, the contributions from string and domain wall decay should 
simply be added to that from vacuum realignment jist . 

In the radiation dominated era under consideration, the space-time met- 
ric is given by: 

-ds^ = -df + R^{t) dx-dx (2) 

where t is the age of the universe, x are co-moving coordinates and R{t) ~ 
is the scale factor. The axion field a{x) satisfies: 

{D^D^^ + ml{t)) a{x) = [d^, + 3^9* - -^V? + ml{t)^ a{x) = (3) 

where ma{t) — ma{T{t)) is the time-dependent axion mass. We will see 
that the axion mass is unimportant as long as ma{t) <^ 1/t, a condition 
which is satisfied throughout except at the very end when t approaches ti . 
The solution of Eq. (|) is a linear superposition of eigenmodes with definite 
co-moving wavevector k: 

a{x,t) = J d^k a{k,t) e*^-^ (4) 

where the a{k,t) satisfy: 

df + ^^dt + ^+ml{t)y{k,t)^0. (5) 

Eqs. and show that the wavelength A(t) — ^R{t) of each mode is 
stretched by the Hubble expansion. There are two qualitatively different 
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regimes in the evolution of each mode depending on whether its wavelength 
is outside (A(t) > t) or inside {X{t) < t) the horizon. 

For \[t) ^ only the first two terms in Eq. (^) are important and the 
most general solution is: 

a{k, t) = ao(fc) + a-i/2(fc)i"'/' • (6) 

Thus, for wavelengths larger than the horizon, each mode goes to a con- 
stant; the axion field is so-called "frozen by causality". 

For \{t) <ti t, the first three terms in Eq. (||) are important. Let a{k, t) = 
R-i{t)ilj{k,t). Then 

{dl^u\t))^{k,t)^0 (7) 
where ^ 2 

Since Co <^ , this regime is characterized by the adiabatic invariant 
ipQ{k,t)uj{t), where ipo{k,t) is the oscillation amplitude of ijj{k,t). Hence 
the most general solution is: 

a(fc,t) = ;^cos(^y dt'uit')^ . (9) 

The energy density and the number density behave respectively as ^ ~ 

WW) 7 iF(t) and j - ip^ j - indicating that the number of 

axions in each mode is conserved. This is as expected because the expansion 
of the universe is adiabatic for X{t)t <C 1. 

Let us define to be the number density, in physical and fre- 

quency space, of axions with wavelength = 77, for oj > t^^. The axion 
number density in physical space is thus: 

na{t) = ^dcu ^{uj,t) , (10) 

whereas the axion energy density is: 

Pa{t)^ [ d0J^{0J,t)uj . (11) 

J^-l duj 

Under the Hubble expansion axion energies redshift according to uj' — lj , 

and volume elements expand according to AV' — AV ("^j > whereas the 
number of axions is conserved mode by mode. Hence: 
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Moreover, the size of for u; ^ j is determined in order of magnitude 
by the fact that the axion field typically varies by fa from one horizon to 
the next. Thus: 



oj— — (w, rjAoj 



^"./i,VlV^i(Va)=^l|. (13) 



du \t' \t 2 ^ ' 2 ^2 



From Eqs. (|T|) and (|T3|), and i? - V*, we have: 

Eq. ([l^ ) holds until the moment the axion acquires mass during the QCD 
phase transition. The critical time is when ma{t) is of order t^^. Let us 
define ti. 

maih)h = 1 . (15) 

ma{T) was obtained H from a calculation of the effects of QCD instantons 
at high temperature H: 



„.„m.4.10-ev(H^)(!|X)' (16, 



3l2( 

V fa J \ T 

when T is near 1 GeV. The relation between T and t follows as usual from 



\2tJ 3 ^ 3 30 ^ ^ 

where is the effective number of thermal degrees of freedom. Af is 
changing near 1 GeV temperature from a value near 60, valid above the 
quark-hadron phase transition, to a value of order 30 below that transition. 
Using A/" ~ 60, one has 

„a,)=0.7 10»^(^)7H;^X-). (18) 

sec Vsecy \ fa J 

which implies: 

The corresponding temperature is: 

Ti ^ IGeV ^^'M . (20) 

We will make the usual assumption |^ that the changes in the axion mass 
are adiabatic [i.e. ^ "'"(f) ^gf^ <C •ma(i)\ after t\ and that the axion to en- 
tropy ratio is constant from t\ till the present. Various ways in which this 
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assumption may be violated are discussed in the papers of ref. We 
also assume that the axions do not convert to some other light axion-like 
particles as discussed in ref. Q . 

The above analysis neglects the non-linear terms associated with the 
self-couplings of the axion. In general, such non-linear terms may cause the 
higher momentum modes to become populated due to spinodal instabilities 
and parametric resonance. However, in the case of the axion field, such 
effects are negligible pO{ . 

We are now ready to discuss the vacuum realignment contributions to 
the cosmological axion energy density. 



2.1. Zero momentum mode 



This contribution is due to the fact that, at time ti, the axion field averaged 
over distances less than ti has in each horizon volume a random value 
between —nfa and +Trfa, whereas the CP conserving, and minimum energy 
density, value is a = 0. In case 1 the average energy density in this "zero 
momentum mode" at time ii is of order: 

pr-'°iti) - lml{h)f!(a'ih)) ^ ^ml{h)f! (case 1) , (21) 

Here a(ti) = a{ti)/fa is the initial misalignment angle. The average 
(a^(ti)) is of order one because a{ti) is randomly and independently chosen 
in each horizon volume. The corresponding average axion number density 
is: 

nr^°iti) = -^,pr'\h) ^ If- (case 1) . (22) 

Since ma{t) is assumed to change adiabatically after ti, the number of 
axions is conserved after that time. Hence: 



n 



vac.O I 



W^^T^f^) (23) 



2ti\R 



The axions in this population are non-relativistic and contribute each 
to the energy. Their average kinetic energy is at all times much less than 
the ambient temperature. 

In case 2 the axion field has been homogenized by inflation and has 
everywhere the same initial value a{ti) — faCe{ti). Hence 

„vac.Ofj-\ ^ fa f „2 



(t)-2^(^) a'ih) (case 2). (24) 



The cold axion density is suppressed if a{ti) happens to be small. The 
probability that it be suppressed by the factor x is of order ^/x. 
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2.2. Higher momentum modes 



This contribution is due to the fact that the axion field has wiggles about 
its average value inside each horizon volume. The axion number density 
and spectrum associated with these wiggles is given in Eq. ( Ilil ) for case 1. 
Integrating over uj > t^^, we find: 











2 ti 





nr'\t)=nr-^ih)^-j^j ^-^^-^j (easel) (25) 

for the contribution from vacuum realignment associated with higher mo- 
mentum modes. The bulk of these axions are non-relativistic after time 
ti and hence contribute nia each to the energy. Note that n™'^'°(i) and 
n™'^'^(t) have the same order of magnitude. 

In case 2, the field is homogenized by inflation and hence n™°'^ is 
very much suppressed. It does not vanish completely because of quantum- 
mechanical fluctuations in the axion field during the inflationary epoch 
[ pl[ . These fluctuations are important if fa approaches the Planck scale. 
However, for fa ^ 10^^ GeV, the resulting n™'^'^ is very small. 



3. String decay 

In the early universe, the axion is essentially massless from the PQ phase 
transition at temperature of order Va to the QCD phase transition at tem- 
perature of order 1 GeV. During that epoch, axion strings are present as 
topological defects, assuming case 1 as we do henceforth. The following 
model describes the relevant dynamics; 

£ = ia^</.ta^0-^(^t^_„2)2 (26) 

where (j) — <j)i + i(f)2 is a complex scalar field. The axion field a{x) is fa 
times the phase of (l){x). The classical solution 

(/.(f) = vaf{p)e'' (27) 

describes a static, straight global string lying on the z-axis. Here {z,p,6) 
are cylindrical coordinates and f{p) is a function which goes from zero at 
p = to one at p = oo over a distance scale of order 6 = {VXva)~^. S is 
called the core size. f{p) may be determined by solving the non-linear field 
equations derived from Eq. (|2^). The energy per unit length of the global 
string is 

H^2Tr [ pdpi|V0p = nvl H^XvaL) (28) 
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where L is a long-distance cutoff. Eq. (|2^) neglects the energy per unit 
length, of order w^, associated with the core of the string. For a network of 
strings with random directions, as would be present in the early universe, 
the infra-red cutoff L is of order the distance between strings. 

At first the strings are stuck in the plasma and are stretched by the 
Hubble expansion. However with time the plasma becomes more dilute 
and below temperature Q 

the strings move freely. String loops of size smaller than the horizon de- 
cay rapidly into axions. The strings that traverse the horizon, called long 
strings, intersect each other with reconnection thus producing sub- horizon 
sized loops which decay efficiently. As a result there is approximately one 
long string per horizon from temperature T, to temperature Ti ~ 1 GeV 
when the axion acquires mass. 

One long axion string per horizon implies the energy density: 

9 

Pstr{t) ^ -^\n{tVa) . (30) 

We are interested in the number density n^*'' it) of axions produced in the 
decay of axion strings because, as we will soon see, these axions become 
non-relativistic soon after time ti and hence contribute each TOq to the 
energy. The equations governing n^*^{t) are: 

dpstr ^rj dpstr^a ^ 

__ ^ _2i/,.,. - (31) 

^ -3i/nf + ^Pflp^ (32) 
dt ° Lu{t) dt ^ ' 

where Loify is defined by: 

1 1 /" dhjdpstr^a ^23) 



uj{t) ^^^f^ J dtduj 



In Eqs. (H-H), '^^7^— is the rate at which energy density gets converted 



dt 

from strings into axions at time i, whereas is the spectrum of axions 

thus emitted. The term —2Hpstr = +Hpstr — '■^Hpstr in Eq. (^ij) takes 
account of the fact that the Hubble expansion both stretches {+Hpstr) 
and dilutes {—3Hpstr) long strings. To obtain nf'^{t) from Eqs. ( ^0[]3l| ), 
we must know what uj(t) characterizes the spectrum dt^dZ"^ '^^ axions 
radiated by cosmic axion strings at time t. 

The axions that are radiated at time t are emitted by strings which 
are bent over a distance scale of order t and are relaxing to lower energy 
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configurations. The main source is closed loops of size t. Two views have 
been put forth concerning the motion and the radiation spectrum of such 
loops. One view [Q, |l^ is that such a loop oscillates many times, with period 
of order t, before it has released its excess energy and that the spectrum 
of radiated axions is concentrated near Let us call this scenario A. A 
second view |||, ^ is that the loop releases its excess energy very quickly 
and that the spectrum of radiated axions ~ ^ with a high energy 

cutoff of order Va and a low energy cutoff of order Let us call this 
scenario B. In scenario A one has uj{t) ~ 2iL and hence, for t < ti: 

nr^W-yln(t'ai) , (34) 
whereas in B one has ~ ^ ln{vat) and hence: 

nfr.^it)^^. (35) 

We have carried out computer simulations |l^ of the motion and decay 
of axion strings with the purpose of deciding between the two possibilities. 
We define the quantity: 

Na.^Jdk'^l (36) 

and compute it as a function of time while the strings decay into axions. 
We have 

_2 ("37) 

where r = iV^"^'/A''™''''' is the factor by which Nax increases during the 
decay of a string loop. Integrating Eqs. ( ^T]) and (|3^), one finds: 
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The string decay contribution is of order r times the contribution from 
misalignment. In scenario A, r is of the order ln(ii/5) ~ 70, whereas in 
scenario B, r is of order one. 

We have performed simulations of two loop geometries: circular loops 
initially at rest, and noncircular loops with angular momentum. The initial 
configurations are set up on large (~ lO*" points) Cartesian grids, and 
then time-evolved using the finite difference equations which follow from 
Eq. (p6|). FFT spectrum analysis of the kinetic and gradient energies during 
the collapse yields Nax{t)- 
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Figure 1. Core position of collapsing loop versus time for A = 0.001 (dotted 
line) and A = 0.004 (solid line). The lattice size is 4096 x 4096. 

3.1. Computer simulations of circular loops 

Because of azimuthal symmetry, circular loops can be studied 'm r — z 
space. By mirror symmetry, the problem can be further reduced to one 
quarter-plane. The static axion field far from the string core is 

a(r,z) = ^l](r,z) (39) 

in the infinite volume hmit, where 17 is the solid angle subtended by the 
loop. We use as initial configuration the outcome of a relaxation routine 
starting with Eq. ( |39| ) outside the core and ^(p) w 0.58(/9/J) exp(i0) within 
the core. Here, p is the distance to the string, and 9 is the winding angle. 
The relaxation and the subsequent dynamical evolution are done with re- 
flective boundary conditions. A step size dt — 0.2 was used for the time 
evolution and the total energy was conserved to better than 1 %. In general 
the loops collapsed at nearly the speed of light without a rebound. For a 
small range of parameters, 80 < Rq/5 < 190, where Rq is the initial loop 
radius, we noticed a small bounce as show in Figure 1. There is a substan- 
tial Lorentz contraction of the string core as it collapses (see Figure 2). A 
lattice effect became evident when the reduced core size becomes compa- 
rable to the lattice spacing. This consists of a "scraping" of the string core 
on the underlying grid, during which the kinetic energy of the string gets 
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Figure 2. Intensity plot of potential energy (contours represent constant poten- 
tial energy) in vicinity of string core for Rq = 2400, A = 0.001 at t = 2600. The 
Lorentz factor is about 4. The arrows represent the axion field. 

Jul 13 1998 Page 1 of 1 

dissipated into high frequency axion radiation. We always choose A small 
enough to avoid this phenomenon. 

Spectrum analysis of the fields was performed by expanding the gradient 
and kinetic energies as 

"^z^ = '^CmnMkmr) sin(A;„2;) 

mn 

Vr(^ = ^CmnJl{kmr) COs(fc„z) 

mn 

<t> = ^CmnJo{kmr) COs{knZ) (40) 

mn 

with the boundary conditions Ji(fcm-Rmax) = and sin(fc„Zniax) = 0. The 
dispersion relationship for axions is given by Wmn = + fc^- Figure 3 

shows the power spcctnim dE /dink displayed in Infc bins of width A In fc = 
0.5, att = and after the collapse at i = 3000. At both times, the spectrum 
exhibits an almost flat plateau, consistent with dE/dk oc 1/fc. The high 
frequency cutoff of the spectrum is increased however after the collapse 
and is associated with the Lorentz contraction of the core. The evolution 
of Nax = J2mni-^mn/kmn) during the loop coUapsc was studied for various 
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Figure 3. Energy spectrum of collapsing loop for _Ro = 2400 and A = 0.001. 
The white (black) histogram represents the spectrum at t = 0{t — 3000). The 
increased high frequency cutoff of the final spectrum is due to the Lorentz con- 
traction of the core. 

values of Rq/S. We observe a marked decrease of Nax by ~ 20% during 
the collapse roughly independent of Rq/S. 



3.2. Computer simulations of rotating loops 

There exists a family of nonintcrsecting ( "Kibble- Turok" ) loops |2^, , 
which have been studied in the context of gauge strings. They are solu- 
tions of the Nambu-Goto equations of motion and have non-zero angular 
momentum. Intercommuting (self-intersection with reconnection) causes 
the loop sizes to shrink, and hence the average energy of radiated axions 
to increase and hence Nax to decrease. Intercommuting favors scenario B 
for these reasons. We picked the Kibble- Turok configuration as an initial 
condition to avoid intercommuting as much as possible and give scenario A 
the best possible chance to get realized. A common loop parametrization 
is given by [ p3| 

-R / 1 
x[a,t) = — I (1 — a) smcT- -I- -a sm3CT_ -f smcr_|_ 
^ \ o 

R / 1 
y{(T,t) — — I — (1 — a) cos (T_ — -a sin 3ct_ — cos ■(/i cos i7_|_ 
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Figure 4. Nax of a circular loop as a function of time for Rq — 2400, and 
A = 0.004 (solid line), A = 0.001 (dashed line), and A = 0.00025 (dotted line). 

z{(j,t) — — y—2y^a{l — a) cosCT- — sini/) coscr+j (41) 

where a± = {a ± t)/R, and a G (0, 2nR) is the length along the loop. For 
a significant subset of the free parameters a € (0, 1), ijj G (— tt, tt), the loop 
never self-intersects. A noteworthy feature in the case of gauge strings is 
the periodic appearance of cusps, where the string velocity momentarily 
reaches the speed of light. 

We performed numerous simulations of rotating loops on a 3D (256'^) 
lattice with periodic boundary conditions. Standard Fourier techniques 
were used for the spectrum analysis, and Nax was computed as a function 
of time using the dispersion relationship 



U!mnp = Y 2(3 - cos km - COS kn - COS fcp) (42) 

Figure 5 shows Nax for various core sizes and constant a, ip. The behavior is 
very similar to that of a non-rotating circular loop with a reduction of Nax 
of about 25%. Figure 6 depicts the energy of the collapsing loop. Clearly, 
the total energy is well conserved. A few percent of the loop energy is 
dissipated as massive radiation, shown here as £^pot- 
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Figure 5. Nax of non- intersecting (" Kibble- Turok") loops as a function of time 
for a = 0.01, = 0, and A = 0.2 (solid line), A = 0.1 (dashed line), and A = 0.0625 
(dotted line). The lattice size is 256^ and R = 72. 

3.3. Conclusions 

In our 3D simulations of rotating loops, r ~ 0.75 for \n{L/5) ~ 4. We take 
L — 2R. In our 2D simulations of circular loops, r ~ 0.8 for ln{L/6) ~ 6. 
The earlier 3D simulations of circular loops by two of us found r ~ 0.8 for 
\n{L/5) ~ 4. We conclude that r is of order one when \n{L/S) ~ 6, and that 
r does not appear to change when ln(L/(5) is increased. Hence we find the 
string decay contribution to the axion cosmological energy density Eq. (js^)] 
to be of the same order of magnitude as that from vacuum realignment 



4. Wall decay 

During the QCD phase transition, when the axion acquires mass, each 
axion string becomes the boundary of N domain walls. If > 1 there 
is a domain wall problem as mentioned in the introduction. We assume 
= 1 here. Each string is therefore attached to one domain wall. The 
domain walls bounded by string are like pancakes or long surfaces with 
holes. A simple statistical argument [|ll| shows that domain walls which 



[Eqs. (Umi)]. 
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Figure 6. Energy of non-intersecting loop. Shown are total, gradient, kinetic, 
and potential energy as a function of time for a = 0.01, tj) = Q, and A = 0.0625. 



are not bounded by string, and therefore close onto themselves like spheres 
or donuts, are very rare. 

When the axion acquires mass, the model of Eq. (^ ) gets modified as 
follows: 

C = ia^</>i5^</'i + \d,hd''h -\{4>l+<t^l- vlf + Tiva<t>i (43) 

The axion mass nia = ^Jr\ for Va 3> ma- We may set 0(a;) — Vae-~°'^^^ 
when restricting ourselves to low energy configurations. The corresponding 
effective Lagrangian is : 





- 1 


cos 







(44) 



Its equation of motion has the well-known Sine-Gordon soliton solution 

a(y) 



27rn -|- 4 tan ^ exp(TOay) 



(45) 



where y is the coordinate perpendicular to the wall and n is any integer. 
Eq. (45) describes a domain wall which interpolates between neighboring 
vacua in the low energy effective theory (^). In the original theory (|4^), 
the domain wall interpolates between the unique vacuum back to that same 
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vacuum going around the bottom of the Mexican hat potential V{<j)^(f>) = 

|((/)t^- 1,2)2 o^^g^ 

The thickness of the wall is of order m^^. The wall energy per unit 
surface is cr = Sniav'^ in the toy model of Eq. (^). In reality the structure 
of axion domain walls is more complicated than in the toy model, because 
not only the axion field but also the neutral pion field is spatially varying 
inside the wall ||2^. When this is taken into account, the (zero temperature) 
wall energy/surface is found to be: 

a ~ 4.2 Um^fa ^ 9 m^/^ (46) 

with fa = Va/N . For N — 1, fa and Va are the same. However, written in 
terms of fa, Eq. ( ^ ) is valid for N ^ 1 also. 

Walls bounded by strings are of course topologically unstable. In ref . [|ll| 
we discussed three decay mechanisms - emission of gravitational waves, fric- 
tion on the surrounding plasma and emission of axions - and concluded that 
decay into non-relativistic axions is very likely the dominant mechanism. 
This is consistent with the computer simulations described below. In the 
simulations the walls decay immediately, i.e. in a time of order their size 
divided by the speed of light, and the average energy of the radiated axions 
is (uJa) 7ma. 

Let be the time when the decay effectively takes place in the early 
universe and let 7 = , be the average Lorentz 7 factor then of the 
axions produced. The density of walls at time ti was estimated ||ll|] to be 
of order 0.7 per horizon volume. Hence the average energy density in walls 
between ti and ^3 is of order 

Pd.w.(t)~(0.7)(9)m,(t):|(^^)'. (47) 

We used Eq. ( ^6| ) and assumed that the energy in walls simply scales as 
niait). After time fs, the number density of axions produced in the decay 
of walls bounded by strings is of order 

Note that the dependence on drops out of our estimate of n'^'^' . 
4-1. Computer simulations 

We have carried out an extensive program of 2D numerical simulations of 



domain walls bounded by strings using the Lagrangian of Eq. (43) in finite 
difference form. We report our results in units where Va — 1 and the lattice 
constant is the unit of length. Thus the wall thickness is 1/ma ~ l/y^ ^^.d 
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the core size is (5 = 1/\/A. In the continuum limit, the dynamics depends 
upon a single critical parameter, niaS = nia/^/X- Large two-dimensional 
grids (~ 4000 x 4000) were initialized with a straight domain wall ini- 
tially at rest or with angular momentum. The initial domain wall was ob- 
tained by overrelaxation starting with the Sine-Gordon ansatz cjji + i(f)2 = 
exp(z tan~^exp(m(jj/)) inside a strip of length D between the string and 
anti-string and the true vacuum (0i = 1, 02 = 0) outside. The string and 
anti-string core were approximated by (t>i+i4>2 = — tanh(.583 r/5) exp(=pi0) 
where (r, 9) are polar coordinates about the core center, and held fixed 
during relaxation. Stable domain walls were obtained for l/{Sma) > 3. A 
first-order in time and second order in space algorithm was used for evo- 
lution with a time step dt — 0.2. The boundary conditions were periodic 
and the total energy was conserved to better than 1%. If the angular mo- 
mentum was nonzero, the initial time derivative (f) was obtained by a finite 
difference over a small time step. 

The evolution of the domain wall was studied for various values of 
^/X/nia, the initial wall length D and the initial velocity v of the strings in 
the direction transverse to the wall. Fig. |^ shows the longitudinal velocity 
of the core as a function of time for the case ml ^ = 400, V\/ma = 10 and 
?; = 0. An important feature is the Lorentz contraction of the core. For 




Figure 7. Speed of the string core as a function of time for the case 1/ma ~ 400, 
y/X/rUa — 10, and v — 0. 
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Figure 8. Position of the string or anti-string core as a function of time for 
l/nia = 1000, A = 0.0002, D = 2896, and v = 0. The string and anti-string 
cores have opposite z. They go through each other and oscillate with decreasing 
amplitude. 

reduced core sizes S/js ^5, where 7^ is the Lorentz factor associated with 
the speed of the string core, there is 'scraping' of the core on the lattice 
accompanied by emission of spurious high frequency radiation. This arti- 
ficial friction eventually balances the wall tension and leads to a terminal 
velocity. In our simulations we always ensured being in the continuum 
limit. 

For a domain wall without rotation {v = 0), the string cores strike 
head-on and go through each other. Several oscillations of decreasing mag- 
nitude generally occur before annihilation. For 7^ ~ 1 the string and 
anti-string coalesce and annihilate one another. For 7., > 2, the strings 
go through each other and regenerate another wall of reduced length. 
The relative oscillation amplitude decreases with decreasing collision ve- 
locity. Fig. H shows the core position as a function of time for the case 
= 1000, VX/rua = 14.3, D = 2896. 

We also investigated the more generic case of a domain wall with angular 
momentum. The strings are similarly accelerated by the wall but string 
and anti-string cores miss each other. The field is displayed in Fig. ^ at 
six time steps for the case — 100, A = 0.01, D — 524 and v — 0.6. No 
oscillation is observed. All energy is converted into axion radiation during 
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a single collapse. 

We performed spectrum analysis of the energy stored in the (t) field as a 
function of time using standard Fourier techniques. The two-dimensional 
Fourier transform is defined by 



fin) 



(49) 
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Figure 10. Energy spectrum at successive time intervals At — 1/ma, with 
Wmax = V8 + m'i and Iw = 18 ln(w/ma)/ ln(TOmax/ma) +2, for the case = 
500, A = 0.0032, D = 2096 and v = 0.25. 



where V = L^Ly. The dispersion law is: 




27rp^ 2TTpy\ 
cos COS I + m^. (50) 



Ly 



Fig. shows the power spectrum dE/dlnuj of the cj) field at various times 
during the decay of a rotating domain wall for the case m^^ — 500, A = 
0.0032, D — 2096, and v — 0.25. Initially, the spectrum is dominated 
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Figure 11. {co) as a function of time for l/rua = 500, D = 2096, v = 0.25 
and A = 0.0004 (solid), 0.0008 (long dash), 0.0016 (short dash), 0.0032 (dot) 
and 0.0064 (dot dash). After the wall has decayed, (uj) is the average energy of 
radiated axions. 



by small wavevectors, k ^ ma- Such a spectrum is characteristic of the 
domain wall. As the wall accelerates the string, the spectrum hardens 
until it becomes roughly 1/k with a long wavelength cutoff of order the 
wall thickness l/nia and a short wavelength cutoff of order the reduced 
core size S/^^. Such a spectrum is characteristic of the moving string. 

Fig. ^ shows the time evolution of {cj)/ma for = 500, u = 

0.25, £> = 2096 and various values of A. By definition, 

(-) = E^p/E? (51) 

P P 

where is the gradient and kinetic energy stored in mode p of the field (p. 
After the domain wall has decayed into axions, (uj) = {uJa}- Fig. |l^ shows 
the time evolution of {Lu)/ma for = 500, A = 0.0016,13 = 2096 and 
various values of v. 

When the angular momentum is low and the core size is big, the strings 
have one or more oscillations. In that case, {io^jraa — 4, which is consis- 
tent with Ref. ||l^ . We believe this regime to be less relevant to wall decay 
in the early universe because it seems unlikely that the angular momen- 
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turn of a wall at the QCD epoch could be small enough for the strings to 
oscillate. 

In the more generic case when no oscillations occur, {uja)/'nna — 7. 
Moreover, we find that {uja)/ma depends upon the critical parameter 
^/X/nia, increasing approximately as the logarithm of that quantity; see 
Fig. |ll|. This is consistent with the time evolution of the energy spectrum, 
described in Fig. 10. For a domain wall, (lo) ^ ma whereas for a moving 
string (w) ~ maln{-</Xva"ts/iTia)- Since we find the decay to proceed in 
two steps: 1) the wall energy is converted into string kinetic energy, and 
2) the strings annihilate without qualitative change in the spectrum, the 
average energy of radiated axions is (uja) ~ mahi{^/Xva/ma)- Assuming 
this is a correct description of the decay process for \/Xva/ma ^ 10^^, then 
{u!a)/ma ^ 60 in axion models of interest. 



4-2. Conclusions 

If we use the value 7 = 7 observed in our computer simulations, Eq. (^) 
becomes ^ 

n^-^ W-f (§) • (52) 
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In that case the contributions from vacuum reahgnment [Eq. ( p3| ) and 
Eq. (^5|)] and domain wall decay are of the same order of magnitude. How- 
ever, we find that 7 rises approximately linearly with ln(VAt;a/ma) over 
the range of \n{^/\va/'ma) investigated in our numerical simulations. If 
this behaviour is extrapolated all the way to \n.{y/\va/rna) — 60, which 
is the value in axion models of interest, then 7 ~ 60. In that case the 
contribution from wall decay is subdominant relative to that from vacuum 
realignment. 



5. The cold axion cosmological energy density 

For case 1, adding the three contributions, we estimate the present cosmo- 
logical energy density in cold axions to be: 

Paito) - 3^ ma {case 1) . (53) 



ti \Ra 

Following 1^ , we may determine the ratio of scale factors Ri /Rq by assum- 
ing the conservation of entropy from time ti till the present. The number of 
effective thermal degrees of freedom at time ti is approximately A/i ~ 61.75. 
We do not include axions in this number because they are decoupled by 
then. Let ^4 be a time (say T4 = A MeV) after the pions and muons an- 
nihilated but before neutrinos decoupled. The number of effective thermal 
degrees of freedom at time ^4 is Af4 = 10.75 with electrons, photons and 
three species of neutrinos contributing. Conservation of entropy implies 
J^iT^Rf = NiTlR\. The neutrinos decouple before e+e" annihilation. 
Therefore, as is well known, the present temperature T^q — 2.735 ''K of 
the cosmic microwave background is related to by: ^r|i?| = 2 T^^ Rl- 
Putting everything together we have 



|l) =0.063(^1 . (54, 



Combining Eqs. (19), (|3) and ( 



Pa{h)-lQ"'''^[-J^S" {easel). (55) 



cm3 V IQi^CeV, 
Dividing by the critical density pc — we find: 

fa 



where h is defined as usual by Hq — h lOOkm/s-Mpc. 
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For case 2, we have 

where a{ti) is the unknown misaUgnment angle. 
6. Pop. I and Pop. II axions 

On the basis of the discussion in sections 2 to 4 we distinguish two kinds 
of cold axions: 

I) axions which were produced by vacuum realignment or string decay 
and which were not hit by moving domain walls. They have typical mo- 
mentum {pi{ti)) ^ aX time ti because they are associated with axion 
field configurations which are inhomogeneous on the horizon scale at that 
time. Their velocity dispersion is of order: 

mah \ R J \ rua J R 

We call these axions population I. 

II) axions which were produced in the decay of domain walls. They 
have typical momentum {pii{t^)) ~ T^TiaC^s) at time when the walls 
effectively decay. Their velocity dispersion is of order: 

Pu{t) - - 10- { (}^\ ^'^ ^ . (59) 

ma R \ nia RiJ \ ma J R 

The factor q = 7 "^°^^*^^ parametrizes our ignorance of the wall decay 
process. We expect q to be of order one but with very large uncertainties. 

There is however a lower bound [0: q > ^ °^ ^ . Since our com- 
puter simulations suggest 7 is in the range 7 to 60, the axions of the second 
population (pop. II) have much larger velocity dispersion than the pop. I 
axions. 

Note that there are axions which were produced by vacuum realignment 
or string decay but were hit by relativistically moving walls at some time 
between ti and ^3. These axions are relativistic just after getting hit and 
therefore are part of pop. II rather than pop. I. 

The fact that there are two populations of cold axions, with widely 
differing velocity dispersion, has interesting implications for the forma- 
tion and evolution of axion miniclusters ||2^, |2^ . We showed that the 
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QCD horizon scale density perturbations in pop. II axions get erased by 
free streaming before the time teq of equaUty between radiation and mat- 
ter whereas those in pop. I axions do not. Hence, whereas pop. I axions 
form miniclusters, pop. II axions form an unclustered component of the ax- 
ion cosmological energy density which guarantees that the signal in direct 
searches for axion dark matter is on all the time. 

The primordial velocity dispersion of pop. II axions may some day be 
measured in a cavity- type axion dark matter detector p8[ . If a signal is 
found in such a detector, the energy spectrum of dark matter axions will 
be measured with great resolution. It has been pointed out that there 
are peaks in the spectrum p9[ | because late infall produces distinct flows, 
each with a characteristic local velocity vector. These peaks are broad- 
ened by the primordial velocity dispersion, given in Eqs. ( ^8|) and ( p9|) 
for pop. I and pop. II axions respectively. These equations give the ve- 
locity dispersion in intergalactic space. When the axions fall onto the 
galaxy their density increases by a factor of order 10'^ and hence, by Li- 
ouville's theorem, their velocity dispersion increases by a factor of order 
10. [Note that this increase is not isotropic in velocity space. Typi- 
cally the velocity dispersion is reduced in the direction longitudinal to 
the flow in the rest frame of the galaxy whereas it is increased in the 
two transverse directions.] The energy dispersion measured on Earth is 
AE = ma/3 A/3 where ~ 10^"^ is the flow velocity in the rest frame of the 

Earth. Hence we find AEj ^ 3- 10~^^ I ^'^^^ I nia for pop. I axions and 

AEjj ^ 10^^^ q ^ '^^ ^ f^a for pop. II. The minimum time required 

to measure AE is (AE)~^. This assumes ideal measurements and also that 
all sources of jitter in the signal not due to primordial velocity dispersion 
can be understood. There is little hope of measuring the primordial ve- 
locity dispersion of pop. I axions since (AEj)^^ ^ 10 years ^ ^^m'^^ ^ 

However [AEu)^^ ~ day q^^ ( ) ' ^^"^ hence it is conceivable 

that the primordial velocity dispersion of pop. II axions will be measured. 
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